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ABSTRACT 

We have identified the third known accretion-powered millisecond pulsar, XTE J0929— 314, with the 
Rossi X-Ray Timing Explorer. The source is a faint, high-Galactic-latitude X-ray transient (d > 5 kpc) 
that was in outburst during 2002 April- June. The 185 Hz (5.4 ms) pulsation had a fractional rms 
amplitude of 3-7% and was generally broad and sinusoidal, although occasionally double-peaked. The 
hard X-ray pulses arrived up to 770 fis earlier than the soft X-ray pulses. The pulsar was spinning down 
at an average rate of v = (—9.2 ± 0.4) x 10~ 14 Hz s _1 ; the spin-down torque may arise from magnetic 
coupling to the accretion disk, a magnetohydrodynamic wind, or gravitational radiation from the rapidly 
spinning pulsar. The pulsations were modulated by a 43.6 min ultracompact binary orbit, yielding the 
smallest measured mass function (2.7 x 10- 7 M Q ) of any stellar binary. The binary parameters imply a 
~ 0.01M Q white dwarf donor and a moderately high inclination. We note that all three known accreting 
millisecond pulsars are X-ray transients in very close binaries with extremely low mass transfer rates. 
This is an important clue to the physics governing whether or not persistent millisecond pulsations are 
detected in low-mass X-ray binaries. 

Subject headings: binaries: close — pulsars: individual (XTE J0929— 314) — stars: neutron — stars: 
low-mass, brown dwarfs — X-rays: binaries 



1. INTRODUCTION 

Accretion-powered millisecond pulsars, the presumed 
progenitors of millisecond radio pulsars, have proven sur- 
prisingly elusive for 20 years. The first example, the X-ray 
transient SAX J1808.4-3658 (P spin =401 Hz, P orb =2 hr), 
was identified as a millisecond pulsar only four years ago 
(Wijnands & van der Klis 1998; Chakrabarty & Morgan 
1998). This only deepened the puzzle of why more ex- 
amples are not known, since SAX J1808.4— 3658 is very 
similar to many of the ~ 50 neutron stars in low-mass X- 
ray binaries (LMXBs) that are not known pulsars (Psaltis 
& Chakrabarty 1999). Earlier this year, a second sys- 
tem was detected, the X-ray transient XTE J1751— 305 
(P spin = 435 Hz, P orb =42.4 min; Markwardt et al. 2002). 
We report here on the discovery of a third example, again 
with a very short binary period. 

The faint X-ray transient XTE J0929-314 (I = 260? 1, 
b = 14?2) was discovered by the All Sky Monitor (ASM) 
on the Rossi X-Ray Timing Explorer (RXTE) in 2002 
April (Remillard 2002). A brief scanning observation with 
RXTE detected persistent 185 Hz pulsations (Remillard, 
Swank, & Strohmayer 2002), and further timing revealed 
a circular, 43. 6- min binary orbit (Galloway et al. 2002). 
The high Galactic latitude makes this source ideal for 
multiwavelength study. Variable optical (Greenhill et al. 
2002; Cacella 2002) and radio (Rupen et al. 2002) counter- 
parts were detected at the X-ray source position measured 
with the Chandra X-Ray Observatory (Juett et al. 2002, 
in preparation), and C III/N III AA4640-4650 and Ha 
A6563 emission lines were reported in the optical spectrum 



(Castro- Tirado et al. 2002). In this Letter, we present a 
detailed analysis of the RXTE observations. 



2. OBSERVATIONS 

XTE J0929— 314 is the faintest new transient discov- 
ered by the RXTE All-Sky Monitor (ASM; Levine et al. 
1996). This 1.5-12 keV instrument consists of three scan- 
ning shadow cameras which provide 90 s exposures of most 
points on the sky every 96 min. XTE J0929— 314 was iden- 
tified using a "deep sky map" technique, in which a map of 
flux residuals is constructed by cross-correlating the pre- 
dicted coded mask pattern against the best-fit data resid- 
uals for each 4' cell in the field of view. Maps are made for 
each ASM camera, using all of the low-background expo- 
sures in a given week. This allows more sensitive searches 
for new sources (as low as ~ 15 mCrab away from the 
Galactic center) compared to the ~ 50 mCrab threshold 
for X-ray error boxes derived from individual camera snap- 
shots. 

After XTE J0929-314 was identified by ASM, we ob- 
tained a series of pointed RXTE observations of the 
source during 2002 May 2 - June 24 (MJD 52396- 
52449). Our analysis is primarily based on data from 
the RXTE Proportional Counter Array (PCA; Jahoda et 
al. 1996), which consists of five identical gas-filled propor- 
tional counter units (PCUs) with a total effective area of 
« 6000 cm 2 and sensitivity to X-ray photons in the 2.5- 
60 keV range. For all our PCA observations, the data were 
collected in GoodXenon mode (in addition to the stan- 
dard data modes), which records the arrival time (1 (is 
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Fig. 1— 2-10 keV ASM and PCA flux of XTE J0929-314 
throughout the 2002 outburst. The 1-d averaged ASM count rate la 
confidence intervals are shown as shaded regions (left y-axis), while 
the PCA integrated flux is overplotted as open squares (right hand 
j/-axcs). Estimated uncertainties on the PCA flux measurements 
are also shown, but are typically smaller than the symbol size. 

resolution) and energy (256 channel resolution) of every 
unrejected photon. We estimated the background flux 
using the "CM" faint-source model for PCA gain epoch 
5 (from 2000 May 13). The photon arrival times at 
the spacecraft were converted to barycentric dynamical 
times (TDB) at the solar system barycenter using the 
Jet Propulsion Laboratory DE-200 solar system ephemcris 
(Standish ct al. 1992) along with a spacecraft ephemcris 
and RXTE fine clock corrections. The position adopted 
was that of the optical counterpart (Greenhill et al. 2002). 
Data from the 20-200 keV High-Energy X-ray Timing Ex- 
periment (HEXTE; Gruber et al. 1996) were also used to 
characterize the hard X-ray spectrum. 

3. ANALYSIS AND RESULTS 

The ASM and PCA intensity history of XTE J0929-314 
is shown in Figure 1. The ASM 1.5-12 keV count 
rate began rising around MJD 52370 and peaked at ss 
2.4 count s" 1 (~31 mCrab) on MJD 52394. The ini- 
tial scanning PCA observation occurred 2 d later, when 
the decay had already begun. Over the following 20 d 
the source faded steadily, holding briefly at 10 mCrab 
on MJD 52420 before declining below the detection limit. 
The last 3er PCA detection was on MJD 52443 at a flux 
of 1.8 x 10~ n ergem^s -1 (2-10 keV). For subsequent 
PCA observations, the flux was below the 3cr detection 
threshold of 7.5 x 10~ 12 ergcm~ 2 s" 1 (2-10 keV; or 1.5 x 
lO" 11 erg cm" 2 s"\ 2-60 keV). The estimated 2-60 keV 
fluence over the entire outburst was 4.2 x 10~ 3 crgcm~ 2 . 
The combined PCA and HEXTE X-ray spectra were con- 
sistent with an absorbed power law+blackbody model, 
with photon index T =1.7-2.0, blackbody temperature 
kT =0.5-0.9 keV, and absorption column density N-r < 
10 22 cm~ 2 . For most of the observations we detected the 
source up to 50 keV. We estimated the bolometric cor- 
rection as the mean ratio of the integrated 2-60 keV to 



2-10 keV flux, 2.34 ± 0.12. No X-ray bursts were detected 
in the data. 

For our timing analysis, we selected photons from the 
top layer of each PCU in the energy range 3-13 keV (ab- 
solute channels 7-30) in order to maximize the pulsed 
signal to noise ratio. We binned the arrival times for 
these photons into 2~ n s (ss0.5 ms) samples. An 185 Hz 
pulsed signal was easily detectable in these data, with 
a fractional rms amplitude of between 3 and 7%. The 
pulse profile in the initial PCA observation on MJD 52396 
contained a weak secondary maximum, but it was broad 
and single-peaked in all subsequent observations except 
for the interval MJD 52425-52430, when a second peak 
was again present (Figure 2). The last 3<r detection of 
pulsations was on MJD 52442, when the source flux was 
7.1 x 10~ n erg cm" 2 s- 1 (2-10 keV). A ~1 Hz quasi- 
periodic oscillation (FWHM=0.21 Hz) with ~5% rms am- 
plitude was detected during most of the observations. 

We investigated the energy-dependence of the pulse 
properties by extracting pulse profiles in nine distinct en- 
ergy bands from one of the longer PCA observations on 
MJD 52408. The rms strength of the pulsation decreased 
from around 5% at 2-4 keV to 2.5% at 10 keV. Com- 
parison of the pulse phase in each energy band indicated 
that significant soft lags were present. For each energy 
band between 2.5 and 10.5 keV the pulse arrival time was 
earlier than in the adjacent lower energy band, with a cu- 
mulative difference of 770 /is. The energy dependence of 
both the rms strength and the phase lag of the pulsation 
appeared to flatten above 10 keV; there were insufficient 
counts to explore the dependence above 20 keV. Similar 
phase lags were observed in the millisecond pulsations from 
SAX J1808.4— 3658, but that source did not show a similar 
energy dependence for the pulsed amplitude (Cui, Morgan, 
& Titarchuk 1998; Ford 2000). We will present a more de- 
tailed analysis of the pulse shape variations elsewhere. 
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Fig. 2.— Pulse profiles of XTE J0929-314 measured on 2002 May 
2, 14 and 31 (MJD 52396, 52408 and 52425). The profile is plotted 
as a solid histogram, with error bars showing the la uncertainties. 
The first two Fourier components are overplotted as dotted lines 
(first harmonic) and dashed lines (second harmonic). Note that the 
phase alignment between the panels is arbitrary. 
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Table 1 

Orbit and Spin Parameters for XTE J0929-3 14 



Parameter 


Value 


Orbital period, P orb (s) 


2614.746(3) 


Projected semimajor axis, a x sin / (light-ms) 


6.290(9) 


Epoch of 90° mean longitude, T n/2 (MJD/TDB) 


52405.49434(1) 


Eccentricity, e 


< 0.007 (2c) 


Spin frequency, vo (Hz) 


185.105254297(9) 


Epoch of spin frequency, to (MJD/TDB) 


52396.0 


Spin frequency derivative, v (Hz s" 1 ) 


-9.2(4) x 10" 14 


x 2 /dof 


1506/907 



The 185 Hz pulsation showed a clear Doppler modula- 
tion, and we used a Fourier frequency history to fit a pre- 
liminary orbital model (Galloway et al. 2002). We then 
applied a more precise analysis using the pulse phase tech- 
niques commonly employed in radio pulsar timing (see, 
e.g., Manchester & Taylor 1977). Adopting our prelimi- 
nary orbital model, we epoch-folded 43.6 min intervals of 
data (thus averaging over periodic phase residuals caused 
by errors in our assumed orbital parameters) and used 
the resulting linear drift in phase residuals to refine our 
knowledge of the pulsar's spin frequency. Then, we com- 
bined this refined spin frequency and the preliminary or- 
bital model to epoch-fold 256 s intervals of data. The 
resulting phase residuals were used to compute differential 
corrections to our Keplerian orbit model (sec, e.g., Deeter, 
Boynton, & Pravdo 1981). No significant eccentricity was 
detectable, but a spin frequency derivative substantially 
improved the fit. Our best-fit orbit and spin parameters 
for the pulsar are given in Table 1. The pulse time delays 
due to the orbit, as well as the effect of fitting a constant- 
frequency model are shown in Figure 3. 



4. DISCUSSION 

From the presence of persistent millisecond pulsations 
over a wide luminosity range in XTE J0929— 314, we can 
infer an upper limit on the pulsar's surface dipole mag- 
netic field strength of B < 1 x 10 9 c?iokpc G from accre- 
tion torque theory (Psaltis & Chakrabarty 1999). This is 
consistent with the expectation that the neutron star is a 
recycled pulsar whose magnetic field has decayed during 
prolonged mass transfer (Bhattacharya & van den Heuvel 
1991). It is interesting to note that the system has bi- 
nary parameters that are extremely similar to those of the 
other recently discovered millisecond X-ray pulsar XTE 
J1751-305 (Markwardt et al. 2002) as well as the slow 
(7.6 s) accreting X-ray pulsar 4U 1626—67 (Middlcditch 
et al. 1981; Schulz et al. 2001), pointing to a similar evo- 
lutionary path. A puzzling aspect is that, unlike the two 
millisecond pulsars, the LMXB 4U 1626—67 has a strong 
(3 x 10 12 G) magnetic field in addition to its slow spin pe- 
riod, which may indicate that its neutron star was formed 
recently through accretion-induced collapse (see Yungel- 
son, Nelemans, & van den Heuvel 2002 and references 
therein). 

We can estimate a crude lower bound on the distance 
to XTE J0929-314 by first noting that the time-averaged 
mass transfer rate driven by gravitational radiation in this 
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Fig. 3.— Pulse timing residuals in XTE J0929-314. Top panel: 
Fit residuals with and without the Keplerian orbit included. Bottom 
panel: Fit residuals for a constant v plus Keplerian orbit. The best 
fit parabola to these residuals is also shown, indicating the need for 
a ;> term in the model. The especially large deviations near MJD 
52425 may be due to pulse shape variations. Error bars represent 
la uncertainties. 



43.6 min binary is (Faulkner 1971) 



Mow = 5 - 5 x 10 



-12 



1.4 Mp 



2/3 



0.01 Mr. 



(1) 

where M x and M c are the neutron star and companion 
masses, respectively. If we assume that the mass accretion 
rate during the outburst is at least as large as Mqw, then 
our non-detection at the end of the outburst requires a dis- 
tance d > 6 kpc. Also, given an outburst recurrence time 
01 ^ 6.5 yr (the time since RXTE was launched), the 2002 
outburst fluence implies a mean mass transfer rate equal 
to Mow for d > 5 kpc. We thus conclude that the source 
lies at least 1.2 kpc above the Galactic plane. Our limits 
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indicate that the mass accretion rate M during the out- 
burst was not more than a few percent of the Eddington 
critical rate. 

The detection of spin-down during the outburst may 
provide an opportunity to test accretion torque theory 
for X-ray pulsars in the low-B regime. The torque on 
a magnetic star from a prograde accretion disk tends to 
spin-up the star for sufficiently high M, while accretion 
will be centrifugally inhibited (the so-called "propeller" 
regime) for sufficiently low M. However, for an interme- 
diate range of M, the positive material torque due to ac- 
cretion may be dominated by a spin-down torque due to 
one of a variety of mechanisms, even while accretion per- 
sists. Two possible spin-down mechanisms include mag- 
netic coupling of the accretion disk and the magnetospherc 
(Ghosh & Lamb 1979, 1991) and expulsion of a centrifu- 
gally driven magnetohydrodynamic wind (Anzer & Borncr 
1980; Arons et al. 1984; Lovelace et al. 1995). In ei- 
ther case, the torque should not exceed the characteris- 
tic accretion torque M\/GM x r co (where r co =50 km is the 
pulsar's corotation radius), which is consistent with our 
measured v for d > 6 kpc. Another possible mechanism 
is gravitational radiation by the rapidly spinning pulsar 
(Bildsten 1998; Andersson, Kokkotas, & Stergioulas 1999; 
Levin 1999). Some of these issues can be explored by ex- 
amining the relation between M and v, which has been 
done extensively for high-i? neutron stars (e.g., Bildsten 
et al. 1997) but never for the low-B case. A detailed 
analysis of the M-v correlation will first require a bet- 
ter understanding of the pulse shape variations, in order 
to limit systematic uncertainties in the frequency history. 
Such work is currently in progress. 

White dwarfs and neutron stars accreting from a 
hydrogen-rich companion cannot evolve to binary periods 
below about 80 min, corresponding to the so-called "period 
minimum" observed for most cataclysmic variables and 
LMXBs (Paczynski & Sienkiewicz 1981; Rappaport, Joss, 
& Webbink 1982). Ultracompact (P orb < 80 min) bina- 
ries like XTE J0929— 314 must therefore have a low-mass, 
hydrogen-depleted (and probably degenerate) donor (Nel- 
son, Rappaport, & Joss 1986). Our measured orbital pa- 
rameters further constrain the nature of the donor in this 



system. The pulsar mass function, fx = 2.7 x 10 -7 M©, is 
the smallest presently known for any stellar binary. This 
mass function gives a minimum companion mass (i — 90°) 
of M c = 0.008 M Q for M x = 1.4 M Q ; it also implies 
M c < 0.03M Q (95% confidence) for a uniform a priori 
distribution in cosi (M x = 2M Q ). The mass-radius rela- 
tion for a Roche-lobe filling donor in a 43.6 min binary 
is R c = 0.04(M C /0.01 MqY^Rq (Faulkner et al. 1972). 
As expected, this has no intersection with the theoretical 
mass-radius relation for very low-mass hydrogen-rich stars 
(i.e. brown dwarfs, Chabrier et al. 2000; cf. Bildsten & 
Chakrabarty 2001). A very low-mass cold helium white 
dwarf (Zapolsky & Salpeter 1969; Nelemans et al. 2001) 
is easily consistent for masses of ~ 0.01 M , implying that 
the system has a relatively high (and thus probable) incli- 
nation. We note that a cold carbon white dwarf is not con- 
sistent with our measured orbital parameters (Lai, Abra- 
hams, & Shapiro 1991). It may be possible for a helium 
dwarf donor to retain a small residual hydrogen content 
(Podsiadlowski, Rappaport, & Pfahl 2002); this is of par- 
ticular interest, given the report of an Ha A6563 emission 
line in the optical spectrum (Castro-Tirado et al. 2002). 

Binary evolution theory predicts that many (if not most) 
of the ~50 known neutron stars in LMXBs are spinning 
at millisecond periods (Bhattacharya & van den Heuvel 
1991). Of these, only three are now known pulsars, with 
persistent millisecond pulsations during their X-ray ac- 
tive states. 1 All three are low-luminosity transients in 
very close binaries, with extremely small time-averaged 
M; evidently, millisecond pulsations are easier to detect 
in such systems. This supports the suggestion that mag- 
netic screening by freshly accreted material may prevent 
the formation of persistent X-ray pulses in NS/LMXBs 
with M above a critical value (Cumming, Zweibel, & Bild- 
sten 2001). 

We are grateful to Jean Swank, Evan Smith, and the 
RXTE operations team at NASA/GSFC for their help in 
scheduling these target-of-opportunity observations. We 
also thank Lars Bildsten, Fred Lamb, Al Levine, Dim- 
itrios Psaltis and Saul Rappaport for useful discussions. 
This work was supported in part by NASA under grant 
NAG 5-9184 and contract NAS 5-30612. 
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